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Context. The y-ray blazar PKS 2155-304 has attracted considerable attention because of its extreme TeV variability characteristics 

during an exceptional flaring period in 2006. Among the observed key findings are (i) a minimum variability timescale as short as 

~ 200 sec and (ii) highly variable TeV emission, which in the frequency interval [10""* Hz, 10"^ Hz] can be described by a log-normal 

distribution and suggests an underlying multiplicative (and not additive) process. 

Aims. Simultaneously accounting for these findings appears difficult within conventional approaches. Following earlier suggestions 

for the TeV blazar Mkn 501 (Rieger & Mannheim 2003), we explore a possible scenario where PKS 2155-304 is supposed to harbor a 

supermassive binary black hole system and where the observed TeV variability is dominated by emission from the less massive black 

hole. 

Methods. We analyze the constraints on the very high energy (VHE) source imposed by the observed variability characteristics and 

the integrated VHE luminosity output, and discuss its implications for a binary black hole system. 

Results. We show that for a secondary mass of men ~ 10'' Mq, fluctuations in the disk accretion rate that feed the jet could account 

for the observed red-noise type variability process down to frequencies of ~ 10"^ Hz. Jet curvature induced by orbital motion, on the 

other hand, could further relax constraints on the intrinsic jet speeds. 

Conclusions. Because a binary system can lead to different (yet not independent) periodicities in different energy bands, a longterm 

(quasi-) periodicity analysis could offer important insights into the real nature of the central engine of PKS 2155-304. 
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1. Introduction 

Active Galactic Nuclei (AGNs) show significant variability over 

a large range of timescales. Usually, a variety of information 

can be extracted from the data with statistical methods based on 

,__! Fourier techniques. In the X-ray domain, for example, variabil- 

J> ity in accreting compact objects has commonly been described 

^\ by means of power spectral densities (PSD), characterizing the 

r — ■ amount of variability power P{a)) as a function of temporal fre- 

OO quencies oj, or timescales l/co (e.g., van der Klis 1997). For 

\l^ (Seyfert) AGNs the resultant PSDs quite often appear to follow 

[^ — . power laws P(oj) oc o)'^, which on long timescales (small fre- 

f^ quencies) are approximately described by j3 close to 1 (flicker 

<^^ noise), but break to a steeper slope (> 2) on timescales shorter 

^~~^ than a break timescale tg. Active Galactic Nuclei thus vary more 

►^ strongly towards lower frequencies (longer timescales). Recent 

• ^H studies have shown that the X-ray PSDs of AGNs can be quali- 

/\^ tatively similar to the high state of black hole X-ray binary sys- 

J-H tems (e.g., McHardy et al. 2004, 2006). Based on these and other 

similarities, AGNs have sometimes been interpreted as scaled-up 

Galactic black hole systems. 

In the very-high-energy (VHE) domain, the experimental situa- 
tion is usually much less favorable. However, one object where 
it became recently possible to employ similar analysis tech- 
niques is the TeV blazar PKS 2155-304 (z = 0.116). Usuafly 
detected only with a low VHE flux of ~ 10% of the Crab neb- 
ula, PKS 2155-304 underwent a dramatic outburst in July 2006, 
with VHE flux levels varying between 1 and 15 Crab units, 
allowing an unprecedented variability analysis (Aharonian et 
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al. 2007; Abramowski et al. 2010; Degrange et al. 2008). The 
Fourier analysis of its VHE light curve of, e.g., MJD 53944 in- 
dicates a red (Brownian) noise-type VHE PSD with an expo- 
nent close to 2 within the frequency range [10 "* Hz, 10"^ Hz] 
(Aharonian et al. 2007). Similar results have been obtained in 
X-rays (1.5-10 keV) with BeppoSAX in 1996 and 1997 (Zhang 
et al. 1999). Moreover, similar to findings in the X-ray domain 
(Zhang et al. 1999, 2005), where linear relations between (abso- 
lute) rms variability amplitude and (mean) X-ray flux have been 
previously reported, a significant rms-flux correlation has been 
observed in the VHE data set (Degrange et al. 2008). This is 
known to be characteristic of a non-linear, log-normal stochas- 
tic process where the relevant, normally distributed variable is 
the logarithm of the flux log(X), and not just the flux itself 
(see Uttley et al. 2005 for a general treatment, and Superina 
et al. 2008 for application to PKS 2155). A log-normal distri- 
bution can be thought of as the result of many multiplicative 
random effects, whereas additive effects would give rise to a 
normal (Gaussian) distribution. Indeed, if the variability would 
be caused merely by an additive stochastic process, the fluxes 
would be normally distributed and no linear rms-flux relation 
would be expected, because all Fourier coefficients would be sta- 
tistically independent. Therefore, the finding that the variations 
on short timescales (determining the rms) decrease in amplitude 
when the long timescale variations (determining the mean flux) 
decrease, indicates that the process driving the variations is a 
multiplicative process (as in a cascade) and not just an additive 
one (e.g., shot-noise), nor one resulting from independent vari- 
ations in many separate regions (Uttley et al. 2005; McHardy 
2008). 
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In accreting galactic black hole (BH) systems these variations 
have been frequently related to small, independent fluctuations 
in the accretion rate, occurring on local viscous timescale fvisc('") 
over a range of disk radii r that are large compared to the inner 
radius of the disk (Lyubarskii 1997; King et al. 2004; Arevalo & 
Uttley 2006). If not damped, these fluctuations can propagate in- 
wards and couple together in a way to produce the multiplicative 
characteristics noted above. Any emission process linked to the 
the innermost region (e.g., as X-ray production or jet launching 
site) may then eventually be modulated over a frequency interval 
ranging from the inverse accretion time near the outer to the one 
at the innermost disk radius, respectively. While this scenario 
sounds very attractive (cf. also Giebels & Degrange 2009 for BL 
Lac), a generalization to a supermassive source like PKS 2155- 
304 seems challenging given the detected VHE minimum vari- 
ability (doubling) timescale of f,, ~ 200 sec. As we show below, 
a consistent approach along this line may require the presence 
of a putative binary BH system in PKS 2155-304 (cf. Rieger & 
Mannheim 2003). It is this conjecture which we wish to explore 
here in more detail. 

2. The origin of variability in PKS 2155-304 

According to the VHE variability analysis for July 2006, there 
are two major findings which need explanation: (i) First, the ex- 
treme short-term variability on a timescale as low as f,, ~ 200 
sec (in the observer's frame^J and (ii) secondly, the log-normal 
distribution of fluxes, which points to an underlying, multiplica- 
tive process. Consider the first finding (i): 
If the luminosity of the host galaxy of PKS 2155-304 is indeed 
as high as reported, i.e., M{R) ^ -24.4 mag (Kotilainen et al. 
1998), its central BH mass is (even if one takes the scatter in the 
Mbh - ibuige relation into account) expected to exceed (McLure 
& Dunlop 2002) 



M, > 2 X 10** Mq 



(1) 



Accounting for the observed rapid variations based on con- 
ventional assumptions has thus proved challenging (see, e.g., 
Levinson 2007 and Begelman et al. 2008 for detailed discus- 
sions): In the galaxy's rest frame, the characteristic scale of any 
disturbance produced by the central engine cannot be smaller 
than its minimum size rg, so that the minimum variability 
timescale for temporal fluctuations of the ejected fluid is rg/c, 
where rg = GM/c^ is the gravitational radius (Levinson 2008). 
Only if parts of the jet flow would collide with a small exter- 
nal disturbance, e.g. a cloud, may smaller variability times get 
imprinted. Yet, in this case the bulk power that can be tapped 
for producing radiation will be correspondingly smaller, and the 
noted PSD shape and lognormal flux distribution may still be in 
need of explanation. If the observed minimum variability would 
indeed reflect the size of the central engine, f,, ~ (1 -H z)rg/c, 
the allowed maximum BH mass would be almost one order of 
magnitude smaller than the one inferred from the host galaxy 
luminosity relation, i.e.. 



M,i4x 10^ ( — - — ) 
1 200 sec/ 



M„ 



(2) 



V 200 sec ^ 

Similar results are obtained if the variability is taken to be pro- 
duced by the collision of shells moving at different speeds (as in 



' As shown by, e.g., Degrange et al. (2008), this value is dependent 
on the telescope's sensitivity, which makes it possible to search for even 
faster variability with future arrays. The currently derived "doubling 
timescale" of ~ 200 sec is therefore best understood as an experimental 
upper limit on the real physical doubling timescale of the source. 



classical internal shock scenarios). In all these cases, the inferred 
minimum source size (in the rest frame of the galaxy) becomes 
independent of any Doppler or Lorentz factorr] Interestingly, 
similar smaU BH masses M, ~ 3x10^ Mq for PKS 2155-304 have 
been previously inferred based on its X-ray variability prop- 
erties (Hayashida et al. 1998; Czerny et al. 2001). Moreover, 
Lachowicz et al. (2009) have recently suggested that the de- 
tection of a possible ~ 4.6 hr X-ray quasi-periodic oscillation 
(QPO) with XMM-Newton in 2006 might be related to a hot spot 
orbiting in the disk close to the innermost stable orbit. If this in- 
terpretation is correct, is would imply a firm upper limit on the 
BH mass of < 3 x 10^ Mq (Schwarzschild) and < 2.1 x IO^Mq 
(extreme Kerr), respectively. 

One obvious possibility to reconcile these apparently divergent 
central mass estimates is related to the putative presence of a 
close binary BH system where the jet that dominates the high- 
energy emission originates from the less massive (secondary) 
BH. In such a case, host galaxy observations may only be indica- 
tive of the total (primary and secondary) BH mass, whereas high- 
energy observations could actually reveal signatures of the less 
massive BH. We note that the possible presence of close binary 
BH systems in blazar-type AGN sources has been repeatedly in- 
voked to account for a variety of observational findings, most 
notably for the detection of mid-and long-term quasi-periodic 
variability (e.g., see Rieger 2007 for a review). 
If the observed VHE output in PKS 2155-305 were be domi- 
nated by emission from the secondary BH, however, its mass 
cannot be supposed to be too small if one also wishes to ac- 
count for the required jet power. Suppose that the maximum 
possible jet power is constrained by the Eddington accretion 
limit Le ~ 10'*^(Mbh/10'^Mq) erg/s. The maximum spin power, 
for example, which can be extracted from the rotational en- 
ergy of a BH (via a Blandford-Znajek-type process) is on the 
order of L, ~ lO'^iMBn/lO^MefiBH/lO^Gf erg/s. With an 
Eddington equipartition magnetic field strength Bh ^ 1/ y/Msn, 
the extractable spin power thus becomes comparable to Le, i.e., 
Ls ~ Le/3. The average VHE gamma-ray flux detected by 
H.E.S.S. implies an isotropic VHE luminosity Liso ~ (5 - 10) x 
10"*'' erg/s (Aharonian et al. 2007). If the radiation is produced 
within two, oppositely directed, conical jets of an opening an- 
gle ~ Q.l rad, the required (minimum) jet power would be 
Lj ~ (1 - 2) X 10^*^(0/0.1)^ erg/s. For a BH accreting at close to 
the Eddington limit (see below) or with maximal spin power, the 
required minimum mass thus becomes 



Ml -10^ Mo (^) , 



(3) 



which would be consistent with the variability constraints (Eq.l2|l 
and be suggestive of a binary mass ratio milnix ^ Msh/M, > 
0.01. 

In reality, the variability constraint of Eq. (J2]i is probably too op- 
timistic, in particular if one wishes to account for the second 
finding (ii) (i.e., the log-normal VHE flux distribution) based 
on fluctuations in the accretion flow that feeds the jet. Consider 
again the fluctuating disk model (Lyubarskii 1997), where fluc- 
tuations of the disk parameters at some radius, which occur on 
local viscous timescale, can lead to variations in the accretion 
rate at smaller radii that are of the flicker- or red-noise type. Then 
the relevant timescale to be employed is the viscous timescale 
fvisc close to the inner radius of the disk, and not the dynamical 



^ Note that the situation could be different if the variations are pro- 
duced within the relativistically moving component itself, e.g., via a 
cascade induced by internal yy-absorption, cf. Aharonian et al. (2008). 
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timescale r^lc. In terms of the a parameter and the disk scale 
height to radius ratio {hji), this timescale can be expressed as 



^vis 



= -(-) 

a\hl 



^'r\"'r. 



(4) 



where a is the viscosity coefficient. Consequently, the inner disk 
needs to be sufficiently thick and viscous ((/i/r)^a ~ 0.1) in order 
to minimize the variability timescale without violating the lumi- 
nosity constraint eq. Q. Similar conditions seem to be required 
to account for the X-ray spectral -timing properties of both AGNs 
and BH X-ray binaries (Arevalo & Uttley 2006). Suppose there- 
fore that in its innermost parts the disk is geometrically thick (as, 
e.g., for an ADAF). Within a more global scenario, one may then 
naturally expect a transition to a standard cooling-dominated 
(geometrically-thin) disk to occur at some radius r^. These tran- 
sition regions are known to be prone to dynamical instabilities 
(e.g., Gracia et al. 2003). Interestingly, on longer timescales the 
X-ray (RXTE) PSD of PKS 2155-304 seems to show a charac- 
teristic roUover-timescale of f^ ~ 1 day (Kataoka et al. 2001). 
If this is indicative of the location of the transition region, i.e., 
fvisc('"tr) ~ fe> ^ reasonable transition radius rti- ~ 50 r^ might be 
inferred (note that for ADAFs (/i/r) ~ 1 throughout). 
In order to account for the obsei-ved VHE variability charac- 
teristics, suppose therefore that fluctuations in the disk accre- 
tion rate on timescales as short as f,, are efficiently transmitted 
to the jet, leading to red noise-type fluctuations in the injec- 
tion rate for Fermi-type particle acceleration. Obviously we will 
only be able to observe emission with red noise-characteristics if 
these signatures do not get blurred by processes occurring on a 
longer timescale within the source: For an observer, flux changes 
will always appear to be convolved and thus dominated by the 
longest timescales (Salvati et al. 1998), i.e.. 



Af, 



obs 



max{f,,, Af,., Afi-ad) 



(5) 



where f,,, Afr and Afraj ^ face are the timescales (as seen by an 
observer) for injection, for photons traveling across the radial 
width of the source, and for the relevant radiative processes, re- 
spectively. For the radial travel time across the source one has 



Atr 



r 
VhC 



(6) 



where F/, > 10 denotes the characteristic bulk Lorentz factor 
and r the radial source dimension in the rest frame of the galaxy 
(note that length contraction only affects the dimension along 
the direction of motion). Typical bulk Lorentz factors inferred 
for the VHE emission region(s) in PKS 2155-305 during the 
2006 high state are in the range of F^ ~ (30 - 50) and per- 
haps even higher (e.g., Begelman et al. 2008). For the model to 
be consistent, one thus requires a radial source dimension that 
satisfies r < 2 X 10'"* (f,/200 sec)(rfc/30) cm. Either the VHE 
region is thus located very close to the BH, or the jet exhibits an 
internal velocity stratification of the (fast) spine - (slow) shear- 
type where the VHE emission is e.g. produced by Compton up- 
scattering within a fast moving spine, similar to the needle-in-jet 
model recently proposed (Ghisellini & Tavecchio 2008). 
Consider next the radiative timescale and assume synchrotron 
losses to dominate the high-energy branch of the energetic parti- 
cle distribution ndY) that is responsible for the up-scattering of 
soft photons to the VHE regime. The radiative timescale thus is 



Afrad - 40 



5x10^ 

r' 



IG 



sec. 



(7) 



Hence, in order to work successfully, the proposed scenario re- 
quires co-moving magnetic field strengths that are sufficiently 
high (B' ~ 1 G) and electrons responsible for up-scattering that 
are sufficiently energetic (Y ~ 10^). These parameters are again 
consistent with recent model fits to the observed spectral energy 
distribution of PKS 2155-304 within the needle-in-jet frame- 
work (Ghisellini & Tavecchio 2008). Obviously, extreme short- 
term variability may for these (and other) reasons not necessarily 
show up at lower energies. As a consequence: If synchrotron ra- 
diation from the fast "needle" component would dominate the 
emission in the optical domain, then the detectable minimum 
variability timescales are expected to be correspondingly longer, 
e.g., on the order of a few to several tens of minutes for the above 
noted parameters. This might be compared with the minimum 
variability timescale of ~ 15 min during an optical obsei-vational 
campaign on PKS 2155-304 in 1995 (Paltani et al. 1997). 

3. Possible implications of a binary BH model 

Following the above analysis, let us suppose that PKS 2155- 
304 harbors a close binary system consisting of a primary BH 
with mass mi ~ 5 x IQ^Mq and a secondary one with mass 
1112 ~ IO^Mq (mass ratio q ~ 0.02). Evidence for a possible 
longterm period of ~ (4 - 7) yr in the optical V band (Fan & 
Lin 2000; cf. also Brinkmann et al. 2000 for soft-X-ray hints) 
might possibly fit well into such a framework: If one requires 
the secondary to be on an orbit that could intersect the accre- 
tion disk around the primary with a typical (maximum) disk size 



of 10^ 



with 



1 and r. 



= 2 Gnii/c^ (Goodman 2003), 
- 2n/Q.t. are constrained to be 



possible Keplerian periods Pj. 

smaller than P^ ^ 45 (1 + qy^^^x^'^(mi /5 x 10** Mq) yr, a condi- 
tion that seems well satisfied for PKS 2155-304. This suggests 
that binary disk interactions could indeed result in optical QPOs 
on the timescale of several years (see Rieger 2007 for more de- 
tails). Accordingly, we may derive an upper limit on the intrinsic 
Keplerian orbital period of the binary 



Pt < 



(l+z) 



r,Opt 

obs 



^ 13 



obs 



yr, 



(8) 



by assuming that the observed optical longterm periodicity is 
caused by the secondary crossing the disk around the primary 
twice per orbital period, a situation expected to occur during 
the initial stage of interaction. In reality, however, this may not 
be fully appropriate because an alignment of the binary's or- 
bital and primary's disk plane typically occurs on comparatively 
short timescales (e.g., Ivanov et al. 1999). One may thus expect 
quasi-coplanar orbits, with the binary possibly surrounded by 
a circumbinary disk, to be much more common. Yet, granted 
that to be the case, periodically modulated mass transfer across 
the gap (see below) could still lead to QPOs in the standard 
disk-dominated infrared-optical regime. Most circumbinary disk 
simulations suggest that the accretion rate is pulsed over one 
orbital period, although occasionally two peaks per period are 
found (e.g., Artymowicz & Lubow 1996; Bate & Bonnell 1997; 
Kley 1999; Hayasaki et al. 2007). This indicates that Eq. ^ still 
provides us with a reasonable constraint on the maximum intrin- 
sic orbital period. If so, then the residual lifetime tg of the sys- 
tem due to emission of gravitational waves (assuming a quasi- 



circular orbit) would be in the range tg oc P 



8/3 



(3xlO''-8xlO^) 



yr. The current separation of the putative binary, on the other 
hand, would be on the order of 

d^6x\Q'H-^\ \ '"' cm- (200-400) r, .(9) 

\13yr/ \5x10'*Mq/ 
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Once the secondary becomes embedded in the outer disk around 
the primary, it starts clearing up an annular gap of a radial size 
Ar ~ {h/d)da^^'^ in the disk as soon as the timescale for gap- 
opening via tidal torques becomes smaller than the timescale 
on which diffusion can refill it (Papaloizou & Lin 1984). It 
was initially thought that this would cut off the binary BHs 
from any mass supply (accretion) of the circumbinary disk. 
However, numerical simulations show that even with significant 
gap, mass can be supplied to the central binary through tidal, 
time-dependent gas streams that penetrate the disk gap, periodi- 
cally approaching and preferentially feeding the secondary (see 
references above). Accordingly, for non-vanishing binary eccen- 
tricity, QPOs with period equal to the orbital one may arise due 
to a time-modulated accretion rate. The simulations indicate in 
particular that a reversal of mass accretion rates can occur (i.e., 
m2 > nil despite 1112 < mi), acting towards equal-mass binaries 
and suggesting that the secondary may become more luminous 
than the primary. If that is the case, then applying variability con- 
straints based on observed emission properties can lead to a spu- 
rious identification of the secondary mass with the (total) central 
BH mass. As suggested above, this may help to explain some of 
the discrepancies in black mass estimates based on high-energy 
emission models and host-galaxy obsei-vations. 
Suppose that the size r^ of the disk around the secondary is com- 
parable to the Roche lobe radius (Eggleton 1983) or the Toomre 
stability radius (whichever is smaller). For 1112 ~ 10^ Mq and 
rj ~ IQ^rg, the viscous timescale at the outer radius fvisc('"rf) ~ 
50 (0.1/a)([r//i]/10)^ yr (standard disk) then most likely ex- 
ceeds the orbital period P^ ~ 10 yr. Any modulation of the mass 
supply to the outer disk at the orbital period may thus be damped 
by the viscosity of the disk, i.e., as the next pulse comes before 
the former was able to propagate in, the amplitude of variability 
becomes small as a result of summation of the pulses (Zdziarski 
et al. 2009). A significant modulation of the jet power (emerging 
from the innermost part of the disk) solely owing to a periodi- 
cally changing mass supply (at the outer disk radius) may there- 
fore not become apparent. On the other hand, if variability in the 
disk occurs on thermal timescale fth ~ o'^iflrgf^^r^jc, optical 
(V band) disk variability on orbital timescale might be possible. 
While a periodically changing mass supply may thus not nec- 
essarily lead to significant periodicities in the jet emission, geo- 
metrical effects due to the orbital motion could introduce them: 
If the jet launched is wrapped by a sufficiently strong magnetic 
field, for example, its overall path can become curved due to the 
orbital motion of the jet-emitting BH (leading to non-ballistic- 
type helical motion). Then the constraints on the required mini- 
mum bulk Lorentz factors of the outflow (Xb,z^ perpendicular to 
the disk plane) might be somewhat relaxed because the effective 
Doppler factor will be time-dependent, D{t), with the strongest 
boosting effects occuiTing whenever the (changing) angle to the 
line of sight becomes smallest (cf. Rieger 2004). This is illus- 
trated in Fig. [T]for a bulk Lorentz factor (in the rest frame of 
the secondary) of 20 and an inclination / (defined as the angle 
between the line of sight and the z-axis, in a frame where the 
orbital plane is in the x-y plane) between 2° and 4°. Obviously, 
considerable changes in the effective Doppler factor could occur, 
which would further influence the detectability of high- and low- 
level (intrinsic) active source stages, i.e., even if the intrinsic flux 
would change only moderately, an order of magnitude differ- 
ence in measured VHE fluxes may become possible. Moreover, 
the motion of components along curved trajectories could possi- 
bly also account for the obsei-ved super-quadratic coiTelation be- 
tween the X-ray and VHE flux (Costamante 2008; Katarzynski 
& Walczewska 2010). Of course, this might only be expected as 




Fig. 1. Time-dependence of the Doppler factor for motion along 
an idealized orbital-driven helical jet path. Chosen parameters 
are: F;,,, = 20, / = 2°, P^. ^ 13 yr (thick solid line), and F^,, = 20, 
/ = 4°, P^ ^ 4 yr (thick dashed line). Also shown is the time- 
dependence of the angle 6{t), the actual angle of the emission 
region to the line of sight (in degree, multiplied by a factor 10) 
for Yb,, = 20, / = 2° and Pj ^ 13 yr. 

long as the lateral width of the jet remains smaller than the sepa- 
ration d of the binary, so that the effects caused by orbital motion 
will not be smeared out. For an expanding jet, this obviously lim- 
its the distance up to which extreme Doppler effects may occur. 
This might partly explain why on larger (radio VLBA) scales, 
only modest Doppler boosting seems to be present (e.g.. Finer & 
Edwards 2004). Clearly, an accurate identification of active VHE 
source stages in PKS 2155-304 could help to further assess the 
influence of differential (geometrical) Doppler effects. 



4. Conclusions 

The observed, extreme VHE variability characteristics of 
PKS 2155-304 provide strong constraints on the physical param- 
eters of its engine. We have suggested that similar to Mkn 501 
(Rieger & Mannheim 2003), the putative presence of a close su- 
permassive binary BH system could allow (i) to reconcile central 
mass estimates based on host galaxy observations (indicative of 
the total primary and secondary mass) with those based on VHE 
y-ray variability (possibly only indicative of the jet-emitting sec- 
ondary), (ii) to account for the observed log-normal variability 
characteristics via accretion disk fluctuations, and (iii) to relax 
constraints on the jet flow velocity. Obviously, an increased in- 
strumental sensitivity in the TeV domain (by, e.g., a CTA-type 
instrument) that may allow us to search for even faster variabil- 
ity and, complementary, an advanced QPO analysis in the optical 
(e.g., Hudec & Basta 2008) could thus be particularly useful to 
further assess the possibility of a binary scenario for PKS 2155- 
304. If coiToborated by further observations, the analogy to X- 
ray binary systems might be closer than initially anticipated. 

Acknowledgements. Discussions with Fehx Aharonian, Amir Levinson and 
Berne Giehels are gratefully acknowledged. 
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